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Abstract 

The decay r~ Itt^tt^Stt^i't- has been studied with the CLEO II detec- 
tor at the Cornell Electron Storage Ring (CESR). The branching fraction 
is measured to be (2.85 ± 0.56 ± 0.51) x 10"^. The result is in good agree- 
ment with the isospin expectation but somewhat below the Conserved- Vector- 
Current (CVC) prediction. We have searched for resonance substructure in 
the decay. Within the statistical precision, the decay is saturated by the 
channels t" TT^2-K^u)i^r, 2t:^ i:^ rjUr ■, and ■k~ 2t:^ 7]^^ . This is the first ob- 
servation of this u decay mode and the branching fraction is measured to be 
(1.89tJ]:^7 ± 0.40) X 10"^- 
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The decay of the r lepton provides a test of the standard model prediction of the hadronic 
weak current. The decay r" 2Ti~TT^'iTi^Vr is related to 37r~27r"'"7r°i/T- and r~ 

TT'^TT^Vr by isospin symmetry and to the cross section for e^e^ Gvr by the Conserved- 
Vector-Current (CVC) hypothesis. The decay can therefore be used to test the isospin 
and CVC predictions. However, these predictions are for the 1 = 1 component of the 
hadronic weak vector current so any J = contribution to the e^e~ cross-section or any 
axial-vector contribution to the r decay through an rj intermediate state must be removed 
before the comparison. The isospin symmetry also relates the relative branching fractions 
among the four possible isospin states 0], 510 (47rp), 330 (3p), 411 (37ru;), and 321 (vrpcij), 
which we denote according to the lowest mass states. The search for possible t] and uj 
substructure will therefore be of particular interest. Any substructure will also be powerful 
in suppressing hadronic background in the measurement of the Vr mass using the six-pion 
decay. Recently, the ALEPH collaboration reported substantial branching fractions for 
the decays r" -> 27r-7r+n7rV^(n > 3) and 27r-7r+37rV^ of (1.1 ± 0.4 ± 0.5) x lO'^ and 
(2.0 ±0.6 ±0.6) X 10~^, respectively. In this Letter, we present a new result for the branching 
fraction of and the first observation of 'k^2'k^ujVt-. We assume that 

all three charged particles in the decay are pions. 

The data used in this analysis have been collected from e^e~ collisions at a center-of- 
mass energy (a/s) of 10.6 GeV with the CLEO II detector at the Cornell Electron Storage 
Ring (CESR). The total integrated luminosity of the data sample is 4.68 fb~^, corresponding 
to the production of 4.27 x 10^ r+r^ events. The CLEO II detector has been described in 
detail elsewhere [Q. 

We select events with four charged tracks and zero net charge. The distance of closest 
approach of each track to the interaction point must be within 0.5 cm transverse to the 
beam and 5 cm along the beam direction. The momentum of each track must be greater 
than 0.02£^beam (-Ebeam = and the polar angle of each track must satisfy | cos6'| < 0.90, 

where 6 is the polar angle with respect to the beam. 

A photon candidate is defined as a crystal cluster with a minimum energy of 40 MeV in 
the barrel region (| cos6'| < 0.80) or 150 MeV in the endcap region (0.80 < | cos6'| < 0.95). 
The cluster must be isolated by at least 30 cm from the projection of any charged track on the 
surface of the calorimeter unless its energy is greater than 250 MeV. In addition, the cluster 
must have a lateral profile of energy deposition consistent with that expected for a photon. 
A subclass of "high-quality" photons is defined to further discriminate against fake photons; 
these photons must have a minimum energy of 150 MeV and pass the isolation requirement 
unless their energy exceeds 250 MeV. All "high-quality" photons must be included in the 
7r° reconstruction. The 7r° candidates are selected based on a requirement on S*^^ = (m^^ — 
mT,o)/a^^, where a^^ is the mass resolution calculated from the energy and angular resolution 
of each photon. We require all vr" candidates to have —3.5 < S^^ < 2.5 and be in the barrel; 
endcap photons are used primarily to veto background events. 

Each event is divided into two hemispheres using the plane perpendicular to the thrust 
axis 1^ , calculated using both charged tracks and photons. There must be one charged track 
in one hemisphere recoiling against three charged tracks in the other (1 vs. 3 topology). In 
the 1-prong hemisphere, the total invariant mass of charged tracks and photons must satisfy 
Ml < 1.0 GeV/c^. We allow up to two high-quality photons in this hemisphere; for the case 
of multiple photons, there must be at least one 7r° candidate. For the 3-prong hemisphere, the 
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FIG. 1. The S^^ distribution for the (a) highest, (b) intermediate, and (c) lowest energy vr'^ (see 
text). The sohd histogram is the sum of the signal Monte Carlo and background (dashed), which 
includes the r migration and hadronic (shaded) background. 

total invariant mass of charged tracks and photons must satisfy M3 < = 1.777 GeV/ c^ [If. 
The magnitude of the total momentum of the particles in the r rest frame, P*, must be less 
than 0.2 GeV/c. In calculating P*, we ignore initial state radiation and assume that the r 
direction is the same as the total momentum vector of the charged tracks and photons in the 
3-prong hemisphere. This requirement selects events with tau-like kinematics, suppressing 
hadronic background and r migration background for which the momentum vector is not a 
good approximation of the r direction. The migration background due to photon conversion 
is further reduced by a cut on the mass of oppositely charged track pairs (Mee), calculated 
assuming the electron mass. In the mass calculation, one of the charged particles must be 
identified as an electron, a particle with a shower energy to momentum ratio in the range, 
0.85 < E/P < 1.10, and, if available, a measured specific ionization loss that is consistent 
with that expected for an electron. Any event with Mge < 120 MeV/c^ is rejected. 

There must be at least six photons in the 3-prong hemisphere, forming three exclusive 7r° 
candidates. If there is more than one combination that satisfies the requirement, the one with 
lowest total = J2^=i[S^'y]i is selected. The S.y.y distributions of the three 7r° candidates, 
classified according to the vr^ energy, are shown in Fig. For each distribution, the S^^ 
of the other two photon pairs must be in their respective signal regions. An enhancement 
at zero is evident in all three distributions, corresponding to the observation of the decay 
27r^7r+37r°z/T-. 

The detection efficiency and background from r migration and hadronic events are 
calculated with a Monte Carlo (MC) technique. We use the KORALB/TAUOLA program [g 
for the r event simulation and the Lund program for hadronic events. The signal decay 
is modeled using a mixture of r~ 2TT~Ti^riVT- 27r~ n^STC^Ur (39%), r~ 7r~27r°?7z/T- — 
27r^7r+37r°z/^ (11%), and 7r^2n^LJi'r 27r^7r+37r°z/T- (50%). The relative mixtures are 
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FIG. 2. The invariant mass spectrum of the 3-prong hemisphere (a) before and (b) after the P* 
cut is imposed. The sohd histogram is the sum of the signal Monte Carlo and background (dashed), 
which includes the r migration and hadronic (shaded) background. 
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FIG. 3. vr+vr vr^ invariant mass spectrum (6 entries/event). The curve shows a fit to the data. 



determined from the measured branching fractions for the first two decays [|T0|. We assume 
that the Sttt] decays proceed through vr/i with a spectral function dominated by the form 

The Sttu; system is modeled using phase space. The 

Since the absolute prediction 



factor of the ai(1260) resonance |TT 



detector response is simulated using the GEANT program |12 



for low-multiplicity hadronic events may be unreliable, the Lund Monte Carlo is used only 
to predict the shape of the 3-prong mass spectrum. The hadronic background is calculated 
by normalizing to the number of data events in the 3-prong hemisphere with M3 > 2.0 GeV 
before the P* cut is imposed. The simulation reproduces the 27r^7r+37r° invariant mass 
spectrum quite well as shown in Fig. |^. An enhancement of events below Mr is evident. The 
signal, background, and detection efficiency are summarized in Table |. Also listed is the 
branching fraction extracted after correcting for the branching fraction of the 1-prong tag of 
(73.0 ±0.3)% §. 

The decay 27r^7r+37r°z/T- can proceed through different intermediate resonances. We 

have recently observed the decays 2Ti~Ti^riVr via 77 — > 77 and 37r° and TT~27r^rii'r 

via 7] —>■ For the latter decay, we also expect an rj signal in the tt+tt^tt^ mass 

spectrum. Figure ^ shows the 7r"'"7r~7r° mass spectrum (6 entries/event), for which the 7r° 
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TABLE I. Summary of signal, background, detection efficiency, and branching fraction. 

Decay Mode 27r~7r"^37r'' ir^ln^uj 

Data 57 19.4+^;^ 

r migration 10.2 ± 2.0 2.3 ± 0.2 

qq background 8.6 ± 2.7 O.OtH 

Eff(%) 2.16 ±0.07 1.65 ±0.09 



-B(xl0-^) 2.85 ±0.56 1^89+^ 



candidates have been kinematically constrained to the nominal 7r° mass. To reduce the 
combinatoric background, we exclude events with a Sn^ mass consistent with that of the t] 
meson. There is an enhancement in the u) mass region corresponding to the first observation 
of the decay 'K^2'K'^ujVr- There is also an indication of a signal for rj tt+tt^tt*^, 

although not statistically significant. To extract the branching fractions, the spectrum is 
fitted using a binned maximum likelihood technique with Gaussians for the uj and rj signals 
and a second-order polynomial background. The uj and rj masses and widths have been 
constrained to the Monte Carlo expectations. The result for r~ — > 7r~27r°a;z/^ is summarized 
in Table |. If we interpret the small r] excess as signal, the branching fraction extracted for 
— s> 'K^2'K^riVr is consistent with the expectation [1^ . 



There are several sources of systematic errors as summarized in Table ||. The uncertainty 
in the photon detection efficiency is estimated by varying the photon selection criteria. The 
systematic error in the decay modeling includes the uncertainties in the branching fractions 
of 271^71^7]^^ and 7r~27r°?7i/,- [|T0| and the modeling of the vr/i and Sttcij systems. The 

latter uncertainty is estimated from the difference in the detection efficiency between the 
default models and the models based on a tt/i phase space and a upoj resonance with mass 
of 1600 MeV/c^ and width of 235 MeV/c^. The uncertainty in the background shape in the 
7r+7r^7r° mass spectrum is estimated by using different orders of polynomial. As a check of 
the hadronic background estimate, the branching fraction for the decay 27r^7r+37r*^z/^ 
has been measured using a lepton tag and the result is consistent with that for the non-lepton 
tag. Including the systematic errors in quadrature, the final results are: 

B{t- 27r"7r+37r°z/^) = (2.85 ± 0.56 ± 0.51) x 10"^ 
B{t- n-2n^uji^r) = (1-89;°;^^ ± 0.40) x 10"^ , 

where the first error is statistical and the second is systematic. 

We can test isospin symmetry [0 by comparing the fractions of the 1 = 1 component of 
the six-pion decay into 27r^7r+37r'^ and 37r^27r+7r° states, 

B{t- 27r-7r+37rV^) 

f27r-7r+37rO " 



B{t- (67r)-i/^) 

and similarly for f3,r-27r+7rO; where B{t^ {67i)^Ut-) is the sum of the three branching 
fractions, B{t~ 27r^7r"'"37r°z/r), B{t~ — >■ 37r~27r"'"7r°^'T-), and B{t~ — >■ ir'biT^i'r). The latter 
branching fraction has not yet been measured. However, we expect the branching fraction 
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TABLE II. Summary of systematic errors (%). 



Mode 


27r-7r+37r° 


7r-27r0a; 


Luminosity 


1 


1 


Cross-section 


1 


1 


T migration 


5 


4 


qq background 


5 


9 


Tracking 


4 


4 


Photon eff. 


15 


15 


Decay model 


4 


8 


Fitting 




5 


Eff (MC Stat.) 


3 


6 


Total 


18 


21 



to be small from CVC using the measured cross section of e~^e Stt+Stt . Sobie 13 and 



Rouge [|1J] have performed an isospin analysis of these decays and concluded that there is 
a discrepancy between the measured decay branching fractions and the isospin expectation. 
However, the authors did not correct for the axial- vector contributions from the decays 

and 7r~27r°?7z/^. We correct for these contributions using our recent measurements 
of the branching fractions, (3.5to:6±0.7) x 10"^ and (1.4±0.6±0.3) x 10"^ [|ig. Figure |shows 



f27r-7r+37r0 VS. f37r-27r+7r0 with our ucw measurement of Byi^T" 27r^7r+37r°z/T-) = (1.41±0.76) x 
10~^ and the world average measurement of i?y(r^ 37r^27r+7r°z/^) = (1.39±0.55) x 10~^, 
where the subscript indicates that this is the vector component of the branching fraction. 
The measurement is in good agreement with the isospin expectation. Because the rj and 
uo intermediate states saturate the decay 27i^7i^37c^h'r, we can also test the isospin 

prediction using our measurement of B{t^ 27r^7r+37r°z/T-) from the decay r~ 'k~2tt^ujVt-. 
With the assumption that B{t~ tt'Stt^z/t-) is negligible, our measurement of B^^ij^ 

(1.68 ± 0.72) X 10"^ is consistent with the isospin expectation as shown 
in Fig. ^ The result prefers the dominance of the [321] (npuj) isospin state over the [411] 
[Siruj) state; we are not able to distinguish the two states from the vrvr^ mass spectrum. Our 
measurement of B{t~ —>■ 27r~7r+37r°z/T-) can also be compared with the CVC predictions by 
Sobie H andEidelman B{t- 27r"7r+37rV^) = (1.9-7.3) x 10"^ and (2.5±0.4) x 10-^ 
respectively. Our result is somewhat below the predictions, indicating that there could be a 
substantial / = contribution from e^e~ rjii^Ti^Ti^ to e^e^ — » 27r+27r^27r'^ that must be 
removed before calculating the CVC predictions. 

In conclusion, the decay r~ — 27r~7r+37r°z/^ has been observed. The branching fraction 
is in a good agreement with the isospin expectation OJT^, but smaller than the ALEPH 
measurement |Q. We have also observed an uj signal in the six-pion system, corresponding 
to the decay n^27T^uJi'r. Within the statistical precision, the contribution from this 

decay together with those |[10| from 27T^n^i]i'r and 7r~27r°r7z/^ saturate the decay 

r~ — > 27r~7r"'"37r°z/^. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 



8 



1160697-010 




0.2 0.4 0.6 0.8 1.0 



FIG. 4. f27r-7r+37r0 VS. fsjr- 27r+7rO • The solid line from the origin uses the CLEO II result on 
Bv{t~ — > 2TT~TT~^37r^h'r) and the world average measurement of Bv{t~ — > 2>'k~2'k^tt^i't) (see text). 
The shaded area indicates the one standard deviation region. The square with diagonal error bar 
uses the CLEO II measurement of B{t~ — > 2-k~ tt^ "iir^ Vr) from the decay r~ — > 'k~ 2it^ ujVr ■ The 
triangular region indicates the isospin expectation. Unitarity requires the fractions to be below the 
dashed line. 

dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A. P. Sloan Foundation. 
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